GRB Hubble diagram and constraints on a A(t)CDM model 
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In previous papers, a cosmological model with constant-rate particle creation and vacuum term 
decaying linearly with the Hubble parameter was shown to lead to a good concordance when tested 
against precise observations: the position of the first peak in the spectrum of anisotropics of the 
cosmic microwave background (CMB), the Hubble diagram for supernovas of type la (SNe la), 
the distribution of large-scale structures (LSS) and the distance to the baryon acoustic oscillations 
(BAO). That model has the same number of parameters as the spatially flat standard model and 
seems to alleviate some observational/theoretical tensions appearing in the later. In this letter 
we complement those tests with 109 gamma ray bursters (GRB), 59 of them with redshifts above 
Cl^ z = 1.4, which permits to extend the Hubble diagram to redshifts up to z ~ 8. For the calibration 

D 1 of the 50 GRBs with z < 1.4 we use the 288 supernovas of the Sloan Digital Sky Survey (SDSS) 

project, calibrated with the MLCS2k2 fitter, less model-dependent than other samples like Union2. 
qq . Our results show a good concordance with the previous tests and, again, less tensions between SNe 

1 la and GRB best fits as compared to the standard model. 
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I. INTRODUCTION 



In the last decade an unprecedented amount of precise observations has led to the consolidation of the so called stan- 
dard cosmological model, with a geometry given by a spatially flat Friedmann-Lemaitre-Robertson- Walker (FLRW) 
metric and an energy content formed by around 5% of baryons, 23% of cold dark matter and 72% associated to a 
cosmological constant A, the component responsible for the present acceleration of the expansion 0. The ACDM 
1 model, together with the standard description of the early hot phase of universe's evolution and the inflationary 
paradigm for the generation of primordial density fluctuations, has been shown robust when tested against the more 
precise observations: i) the spectrum of anisotropies of the CMB temperature 0; ii) the mass power spectrum of LSS 
distribution Q ; iii) the Hubble diagram for supernovas of type la Jj, 0, H| ; and (in a still lower precision level) the 
distance scale to the BAO peak [tJ. The two later are background tests of the redshift-distancc relation predicted 
by the model, and involve relatively low redshifts (up to z ~ 2). The two former also test predictions at the per- 
, turbative level and up to the redshift of last scattering and beyond. For example, the matter density at the time of 
\JT) • matter-radiation equality (around z w 3300 for the ACDM model) is essential to determine the turn-over of the mass 
^sO \ power spectrum and hence its present shape. In the case of the CMB anisotropy spectrum, many free parameters are 
■ available in the data fitting. If we are just interested in fixing background parameters like the present relative matter 
density for example, a simplified test is given by fitting the position of the first acoustic peak in the spectrum. 
CN| ■ In spite of the relative success of the standard model in getting a good concordance between those observations, 
some theoretical and observational problems have stimulated the study of alternative models of dark energy. Among 
the theoretical problems usually mentioned we can cite the absence of a microphysical explanation for the tiny value 
of the cosmological constant and the approximate coincidence between A and the present matter density. On the 
observational level, the celebrated concordance is actually dependent on the sample of SNe la used in the joint analysis 
Q. Among the recent compilations, the Union2 Q is the biggest of the samples we will consider here and, alone, 
predicts for the ACDM a present matter density parameter f2 m o ~ 0.27. However, we should have in mind that the 
Union2 datasct is calibrated with the Salt2 fitter, which makes use of a fiducial ACDM model for including high-z 
supernovas in the calibration. Therefore, that sample is not model-independent and the test should be viewed as 
rather a test of consistence. 

Model-independent SNe samples can be obtained by using instead the MLCS2k2 fitter, which uses only low and 
intermediate redshifts in the calibration H. Although this may lead to other kinds of limitations, it is a must if we 
want to test the robustness of the concordance. By using, for instance, the Constitution and SDSS samples calibrated 
with that fitter, we obtain fi m0 ~ 0.33 and fi m o ~ 0.40, respectively (see e.g. table II of reference [lj]), values in 
clearly tension with that obtained from LSS. For example, the best-fit value obtained from the analysis of the 2dFGRS 
data for the mass power spectrum is il m o ~ 0.23 [Toj . A similar result, fi m o ~ 0.22, is obtained from the analysis of 
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the shape of the SDSS DR5 galaxy power spectrum when only large scales (0.01 < k < 0.06 /iMpc -1 ) are included 
[Til ]. This tension is also clear when we look the reduced x 2 's of the joint analysis including the four tests referred 
above (see Table I below). 

Such a tension may be just the result of systematics in the samples or in the calibration. But it may also have 
some physics behind it. As we have commented above, the mass power spectrum is particularly sensitive to the 
matter density at the time of matter-radiation equality. While the SNe la distance relations are tested for relatively 
low redshifts. Therefore, the mentioned tension may be suggesting a late-time process of particle creation from the 
vacuum, in such a way that we have at early times the matter density needed to fit the LSS observations, but a 
present density above that expected if no matter has been created. 

In order to verify this hypothesis, the data must be fitted by a model with particle creation. Since the radiation 
epoch must remain unaffected, no photons should be produced. In the same way, production of baryons is not allowed 
if we want to maintain the main features of CMB. Therefore, only dark particles will be created in this model. 
Furthermore, since the observed dark matter is cold, the created matter is pressureless. The simpler way of realizing 
this phenomenological recipe is with a constant rate of creation. As we will see, the corresponding model is uniquely 
defined, has the same number of free parameters as the standard model and is not reducible to the later. The data 
analysis with this model leads to a good concordance. With the Union2 sample the fit is as good as in the standard 
case, while the use of the Constitution and SDSS samples leads to a very good agreement with the value derived 
from the LSS analysis in this model, fi m o ~ 0.45. We will also see that the backreaction associated to the matter 
creation implies a time-dependent vacuum density, linearly proportional to the Hubble parameter, in order to respect 
the conservation of total energy. 

In addition to the main observational tests referred above, several others can be done, despite the less precision 
involved. Among them we can cite, for example, weak gravitational lensing, gas mass fraction in clusters, redshift- 
distancc relations for compact x-ray sources, the redshift-age relation for high-z objects, and the peculiar velocity 
field of galaxies. For elucidating the problem we are discussing here, i.e., the tension between early and late-time 
predictions for the matter density, an interesting complementary test is provided by gamma ray bursts (GRB) [l2|, 
since they can extend the Hubble diagram up to redshifts z ks 8 [HI, [l5[. Actually, the use of GRBs as standard 
candles is quite controversial, as we will discuss in the Section III. In this paper we do assume that GRBs can be used 
as cosmological probes, following references [H, EH • O nce more, any strong model dependence should be avoided 
and, for this goal, new GRB samples will be generated in this letter, calibrated with help of the SDSS (MLCS2k2) 
SNe la compilation discussed above. 



II. THE MODEL 



Consider the possibility of creation of non-relativistic dark particles at a constant rate T. The balance equation for 
this process (some times called the Boltzmann equation) is given by 

4,| (A) -ix (i) 

where n is the particle number density at a given time. By multiplying by the created particle mass M, it can be 
rewritten as 

Pm + 3Hp m = Tp m , (2) 

where p m = nM is the matter density. 

In addition, let us take the Friedmann equation for the spatially flat case, 

P = Pm +A = 3H 2 , (3) 

where p m is the matter density and A is the energy density associated to a vacuum term. Combining the two above 
equations and using for the vacuum term the equation of state p\ = —A, we can obtain 

p + 3H(p + p) = (4) 

(where p and p refer to the total density and pressure, respectively), provided we take A = 2TH + Ao, where Ao is a 
constant of integration. Since there is no microphysical scale defining such a constant, we will take Ao = 0. Therefore, 
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TABLE I: Limits to Q ro0 (SNe+ CMB + BAO+LSS) [2g|. 

A(t)CDM ACDM 

Union2 (SALT2) 0.420±g;g?o 1-063 0.235 ± 0.011 1.027 

SDSS (MLCS2k2) 0.450lg;g^ 0.842 0.260±g;gH 1.231 

Constitution (MLCS2k2[17]) 0.450±g;^| 1.057 0.270 ± 0.013 1.384 

"Error bars stand for 2a 

we have 1 

A = 2TH, (5) 

and equation Q can be rewritten as 

p m + 3Hp m = -A. (6) 

We can see that the creation of pressureless matter is only possible with a concomitant time- variation of the vacuum 
density, which acts as a source of particles in the conservation equation, Eq. ([6]), decaying linearly with the Hubble 
parameter H . 

Dividing Eq. (Jjjjl by 3H 2 and using A/(3H 2 ) = 1 - Q m , we obtain 

r = |(i-o ro )/f, (7) 

which is valid for any (late) time. For the present time (with Sl m o ~ 1/3) we have T Hq. On the other hand, in 
the future de Sitter limit (Q TO — s- 0), it follows that T « (3/2)Hds- That is, apart from a constant factor, the creation 
rate is equal to the Gibbons-Hawking temperature associated to the de Sitter horizon fl9j . 
With A = 2TH it is easy to obtain the Hubble function [2(J UK 

J- = i - n TO0 + n m0 (i + zf' 2 , (8) 

where £l m o and Hq (the present values of the relative matter density and of the Hubble parameter, respectively) are, 
like in the standard model, the only free parameters to be adjusted. In the asymptotic limit z — > — 1 it leads to a 
constant H, that is, to a de Sitter universe. On the other hand, for early times (high z), we have H « H n mQ z 3 / 2 , 
which, through the Friedmann equation, leads to p m = 3HQCt 2 n0 z 3 . Therefore, at early times the energy density scales 
with z 3 , that is, with a -3 , as in the standard model. However, in the later we would have p m — 3i?Qf2 m oz 3 . The 
extra factor Q m o (< 1) in the former is related to the late-time matter creation. In order to have the same amount 
of matter today, we should have less matter at higher redshifts. Reversely, if we want the same amount of matter at 
high z (in order to preserve the observed features of CMB and LSS), we expect a present matter density higher than 
in the standard model. 

This is confirmed when we perform precise observational tests with the present model, at both the background and 
perturbative levels. In previous publications [3, [2(| HH, [H, [23|, [U, [25| the model has been tested against observations 
of the Hubble diagram for type la supernovas, the position of the first peak in the CMB spectrum of anisotropies, the 
power spectrum of large scale structures and the distance to baryon acoustic oscillations. The concordance is quite 
good and, depending on the supernovae sample used, it is actually better than in the ACDM case, as shown in Table 
I (with Hq marginalized). 

We can see in the table that, when the Union2 compilation is used in the joint analysis, similar reduced x 2 's are 
obtained for both models. However, Union2 data are calibrated with the Salt2 fitter by making use of a fiducial ACDM 



1 Strictly speaking, this result is only exact if we neglect the conserved baryons in the balance equations. Since baryons represent only 
about 5% of the total energy content, relation (0 can be considered a good approximation. 
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model in the calibration process. A model-independent analysis is possible if we consider instead the MLCS2k2 fitter 
in the calibration, as is the case of the SDSS and Constitution samples in Table I. With these compilations we obtain 
a better concordance with the present model. As anticipated above, the present matter density is higher than in the 
standard case, with O m o ~ 0.45. 



III. THE GRB SAMPLES 



The luminosity of GRBs appears to be correlated with their temporal and spectral properties and, although these 
correlations are not yet fully understood from first principles, their existence has naturally suggested the use of GRBs 
as distance indicators, offering a possible route to probe the expansion history of the Universe up to rcdshifts z ~ 8 

The typical spectrum of the prompt emission of GRBs can be expressed as an exponentially connected broken 
power-law, the so-called Band function [26|. Then we can determine the spectral peak energy E p , the photon energy 
at which the vF v is brightest. 

The two most used correlations are the Amati relation E p — E isoi relating the rest frame energy of the spectra E p 
and the emitted isotropic energy Ei SO , and the Ghirlanda relation E 7 — E Pl relating E p with E 7 = Ei so (l — cos 8 'jet)- 
The latter correlation is the tightest of the GRBs calibration relations. However, to be included in this relation 
the GRB afterglow must have an observed jet break in its light curve, thus only a fraction of the observed events 
can contribute for establishing this relation. Therefore, the Amati relation is more statistically reliable as it can be 
calibrated with a greater number of events. 

We shall apply here Amati's power law expression JUI^, relating E p = E p ^ s ( \ + z) with the isotropic equivalent 
radiated energy Ei SO , defined as 

E lso = 4ird 2 L S bo i {l + z)-\ (9) 

where Sboio is the bolometric fluence of gamma rays in the GRB at a redshift z and g?l is the luminosity distance to 
it. Similarly to [l3j . we will rewrite the Amati relation as 

lo g ^i = A + 61og-^-. (10) 
erg 300keV 

GRBs can be considered as distances indicators if they can be calibrated at low rcdshifts. However, a few data 
are available at low redshifts. To mend this lack of data most calibrations take for granted a particular cosmological 
model. In order to avoid this circularity problem here we follow the method proposed in [2t| , which consists in using 
SNe la as a distance ladder to calibrate the GRBs. Note that with the well-known relation 

one can convert the distance moduli \i into luminosity distances cLl (in units of Mpc). 

Since the distance moduli for the SNe la are known, a cubic interpolation is performed to determine the parameters 
A and b in the Amati relation for the low- redshift GRBs (z < 1.4). Then, the distance moduli for the high-z GRBs 
are obtained and they can be used as standard candles without the circularity problem. It is worth mentioning that 
the calibration of GRBs is still a quite controversial subject. Indeed, even the Amati relation has been contested by 
some authors (see, for instance, [301, where it is argued that the Amati relation could be an artifact of selection effects 
in both the burst population and the detector). With this cautionary remark, we proceed with the calibration. 

We perform the cubic interpolation by two methods: splines and Hermite. A spline is a polynomial between each 
pair of table points, but one whose coefficients are determined "slightly" nonlocally. The nonlocality is designed to 
guarantee global smoothness in the interpolated function up to some order of derivatives. The cubic spline interpo- 
lation produces an interpolated function that is continuous through the second derivative [3l| . On the other side, 
Hermite interpolation is also a method of interpolating data-points as a polynomial function. The generated Hermite 
polynomial is closely related to the Newton polynomial, in that both are derived from the calculation of divided 
differences. Hermite interpolation (first-derivative continuous) requires four points so that it can achieve a higher 
degree of continuity (32[. 

So, by using at low redshifts the 288 SDSS SNe la data (33|, we derived the distance moduli for the 50 low- redshift 
(z < 1.4) GRBs by cubic-splines and cubic-Hermite interpolations. From Eq. © with the corresponding Sboio and 
luminosity distances cIl, we then derived Ei SO for these 50 GRBs. Furthermore, with the corresponding E p we found 
the best fit for the Amati relation (flOl): 



A = 52.8076 ± 0.0683, b = 1.6909 ± 0.1156, 



(12) 
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FIG. 1: Solid lines are the 68.3%, 95.4% and 99.73% C.L. constraints on the free parameters of the A(t)CDM model. The 
black dot is the best fit for each GRB sample. All GRB samples have been calibrated with the Union2 SNe la data set (closed 
elliptical contours). 



for the cubic-splines interpolation, and 

A = 52.8960 ±0.0874, b = 1.7851 ± 0.1480, (13) 

for the cubic- Hermite one. Extrapolating these calibrated Amati relations to the 59 high-redshift (z > 1.4) GRBs 
given in Table 2 of [HI , the distance moduli fj, are obtained. Only these new 59 datapoints are then used in order to 
test the cosmological model 2 . 



IV. RESULTS AND DISCUSSION 



In the last section, assuming the SDSS (MLCS2k2) SNe la sample as our calibration source at low redshifts, we 
have produced two new GRB data sets based on the Hermite and the spline interpolation methods. Before procee ding 
with a statistical analysis with our new samples we consider three previous GRB data sets: i) the NPZ sample 1341] 
with 27 GRBs at z < 1.4 and the remaining 42 objects in the range 1.4 < z < 6.6; ii) the Hymnium sample fla ). 
with 50 low-redshift data points and 59 GRBs in the range 1.4 < z < 8.2; hi) the MB sample |18|, that makes use of 
the same GRB compilation as the Hymnium data, but using a different interpolation for the 50 low-redshift GRBs, 
leading to a different extrapolated Hubble diagram at high z. These samples share the same SNe la calibrating source, 
the Union2 data. A direct comparison of the A(t)CDM model with such samples retains information from the Union2 
data set and, consequently, an imprint of the ACDM model. Nevertheless, for the sake of comparison, we will perform 
a statistical analysis with all the above GRB samples. 

Assuming a set of free parameters (fl m .o, h), where Ho = lOO/i Km s _1 Mpc _1 , the agreement between theory and 
observation is measured by minimizing the x 2 > 

N (n th - n obs ) 2 

x 2 = X 2 (^ m0 , &) = £ [ - ft j , (14) 

i=l U * 

where [i th and fi obs are the theoretical and observed values of the distance moduli, respectively, and a means the 
error for each one of the N data points. The 68.3%, 95.4% and 99.73% confidence levels (C.L.) are determined by 
Ax 2 = X 2 - xliin < 2 - 3 ' 6 - 17 and u - 8 > respectively. 

Fig. [T] shows the constraints on the free parameters (Q m o,h) for the A(t)CDM model using the above mentioned 
NPZ (left panel), Hymnium (central panel) and MB (right panel) samples. The contours for the Union2 SNe la 
sample are also shown. In spite of the large dispersion observed in the GRB samples one can infer from their best fits 



2 Both samples are available under request. 
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FIG. 2: Solid lines are the 68.3%, 95.4% and 99.73% C.L. constraints on the free parameters of the ACDM model. The black 
dot is the best fit for each GRB sample. All GRB samples have been calibrated with the Union2 SNe la data set (closed 
elliptical contours). 




FIG. 3: Constraints on the free parameters of the A(t)CDM (ACDM) model (left (right) panel). The 59 GRBs at 2 > 1.4 
were calibrated with the SDSS (MLCS2k2) SNe la data set using cubic interpolations with hermite and spline methods. The 
GRB contours correspond to the 68.3% C.L. (with best fit at the black dot) for each calibration. The ellipses correspond to 
68.3%, 95.4% and 99.73% C.L. constraints from the SDSS SNe la sample. 



[NPZ: (0.463,0.727), Hymnium: (0.422,0.708), MB: (0.439,0.757)] that they agree with larger values for the amount 
of dark matter and the Hubble constant as compared with the Union2 best fit (O m o = 0.337, h = 0.697). A similar 
analysis for the ACDM model is displayed in Fig. O Note that in this case very large values for the Hubble constant 
are predicted. 

With the new GRB samples of the last section we obtain the results shown in Fig. [3] for the A(t)CDM (left panel) 
and ACDM (right panel) models. The 68.3% C.L. (with best fit at the black dot) are shown for each calibration. Solid 
(dashed) open contours correspond to the cubic Hermite (spline) interpolation method. Again, the large dispersion 
does not allow a definite conclusion, except that both models are in accordance with data within la level. However, a 
look at the best-fit points is suggestive, specially if combined with the better posed results discussed in the Introduction. 
In the ACDM case the matter density parameter is in agreement with that resulted from the LSS analysis and in 
conflict with the larger value obtained with SDSS SNe la alone. If we consider instead the scenario with particle 
creation, we have again a very good concordance between the GRB and SDSS SNe la predictions for the present 
matter density, with both, Hermite and spline methods. If we also take into consideration the results for the Hubble 



7 



constant, the concordance is particularly good with the Hcrmite interpolation. In none of the models the best-fit 
values for the matter density change when h is marginalized. 

In conclusion, we can see that the creation of particles from vacuum, something allowed by and generally expected 
from quantum field theories in curved spacetimes, is not only in agreement with current cosmological observations 
but is indeed favored by them when model-independent datasets are used. In the present model, where pressureless 
dark matter is produced with a constant rate, the vacuum density scales linearly with H in order to conserve the total 
energy. Interestingly enough, this time-dependence for A is precisely that derived when we estimate the renormalized 
density of the QCD vacuum condensate in the FLRW spacetime (see, for instance, 25( and references therein). In 
this case the creation rate is T k to 3 , where to 150 MeV is the energy scale of the QCD chiral phase transition. 
This value for T leads to a remarkable coincidence between A and the observed dark energy density. Furthermore, as 
discussed in [25j , the energy flux from vacuum to the matter sector also aliviates the cosmic coincidence problem. We 
think that all these theoretical and observational positive results are enough reasons to consider that scenario and to 
test it against a wider set of observations. 
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